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Abstract 

Theta vacuum effects on the QCD phase structure in the fi-T plane are studied by using the Polyakov- 
loop extended Nambu-Jona-Lasinio model and its extension, where fi is the quark chemical potential and T 
is temperature, respectively. As the parameter 9 of the theta vacuum increases, the chiral transition becomes 
stronger. For large 9, it eventually becomes first order even at zero fi. 
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I. INTRODUCTION 



Violations of parity (P), charge conjugation (C) and charge-parity symmetries (CP) are im- 
portant subjects in particle and nuclear_physics. For example, the strong CP problem is a long- 
standing puzzle; see for example Ref. [1] for a review of this problem. Lorentz and gauge invari- 
ance allow the Quantum Chromodynamics (QCD) action to have a term 
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of the topological charge, where F^j, is the field strength of gluon. The parameter 9 can take 
any arbitrary value between — tt and tt, where 9 = —n is identical with 9 = n. Nevertheless, 
experiment indicates 16*1 < 3 x 10"^° Iz, sj]. Since 9 is P-odd (CP-odd), P (CP) is then preserved 
for = and ±7r, but explicitly broken for other 9. Why is 9 so small ? This is the so-called 
strong CP problem. 

For zero temperature (T) and zero quark-chemical potential (/x), P is conserved at 6* = 0, as 
Vafa and Witten showed Meanwhile, P is spontaneously broken at 6* = tt, as Dashen [5] 
and Witten pointed out. This is the so-called Dashen phenomena. Since the spontaneous 
P violation is a nonperturbative phenomenon, the phenomenon was so far studied mainly with 



1211 . the Nambu-Jona-Lasinio (NJL) 



the effective model such as the chiral perturbation theory [|7|- 

model il3l- ll5|l and the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model |16|] . 

For T higher than the QCD scale Aqcd, there is a possibility that a finite 9, depending on 
spacetime coordinates (t, x), is effectively induced, since sphalerons are so activated as to jump 
over the potential barrier between the different degenerate ground states nlln . If so, P and CP 
symmetries can be violated locally in high-energy heavy-ion collisions or the early universe at 
T ^ Aqcd- Actually, it is argued in Refs. [|l2L llSQ that 9 may be of order one during the QCD 



phase transition in the early universe, while it vanishes at the present epoch [[19 



2M- This finite 



9 may be a new source of very large CP violation in the Universe and may be a crucial missing 
element for solving the puzzle of baryogenesis. 

Furthermore, this effective 9{t, x) deviates the total number of particles plus antiparticles with 
right-handed helicity from that with left-handed helicity. The magnetic field, formed in the early 
stage of heavy-ion collision, will lift the degeneracy in spin depending on the charge of particle. 
As a consequence of this fact, an electromagnetic current is generated along the magnetic field, 
since particles with right-handed helicity move opposite to antiparticles with right-handed helicity. 
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24l-l26ll. CME may explain the charge 



This is the so-called chiral magnetic effect (CME) [|l8 
separations observed in the recent STAR results Vi21n . Thus, the thermal system with nonzero 9 is 
quite interesting. 

In this letter, we study effects of the theta-vacuum on QCD phase diagram by using the two- 
flavor PNJL model [|l6L I28l-l44ll and its extension; see Ref. nloA and references therein for further 
information on the PNJL model. Particularly, our attention is focused on the two-flavor phase 
diagram in the fi-T plane, where yU is the quark chemical potential. In our previous work Ref. [|l6n . 
as a theoretical interest, we investigated spontaneous P and C violations at finite 9 and imaginary 
/i. As physical phenomena that may occur in the early universe or the high-energy heavy-ion 
collisions, we here examine the chiral and P symmetry restorations at finite T and 9 and real fi. 

This paper is organized as follows. In section II, the PNJL model and its extension are explained 
briefly. In section III, the numerical results are shown. Section IV is devoted to summary. 



II. PNJL MODEL 



Pioneering work on the parity violation and its restoration in the framework of the NIL 



model 113 
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led this issue extensively [|14l . 



45I-I47II was done by Fujihara, Inagaki and Kimura tXM- Boer and Boomsma stud- 



ISO. In Ref. [la], we have extended their formalism based on the 
NIL model to that on the PNJL model. The two-flavor (Nj = 2) PNJL Lagrangian with the 
6'-dependent anomaly term is 



£ = q{i-f^D'' -m)q-l({<P[A],<P[A]*,T) 

3 

+ GiJ2 f + (^^75rag)'] + 8G2 [e'' det {qnqL) + e"'' det {q^qR 



(2) 



a=0 



where q = {u, d) denotes the two-flavor quark field, m does the current quark-mass ma- 
trix diag(m„, m^i), tq is the 2x2 unit matrix, Ta{a = 1,2,3) are the Pauli matrices and 
= d" + iA" - ifi6^. The field A" is defined as A" = S^gA'l^ with the gauge field A^, 
the Gell-Mann matrix A„ and the gauge coupling g. In the NIL sector, Gi denotes the coupling 
constant of the scalar and pseudoscalar-type four-quark interaction^ and G2 is the coupling con- 



stant of the Kobayashi-Maskawa-'t Hooft determinant interaction 1148 . 



49h the matrix indices of 



which run in the flavor space. The Polyakov potential U, defined later in (1191 ), is a function of the 
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Polyakov loop ^ and its Hermitian conjugate <P*, 



<P = -^TtL, (P* = -^TiL^ (3) 



with 



L(x) = V exp 



i I dTAJx,T 







(4) 



where V is the path ordering and A4 = iA^. In the chiral limit (r/iu = ma = 0), the Lagrangian 
density has the exact SU{Nf)L x S'f/(A^/)R x U{1)^ x SU{3)c symmetry. The f/(l)A symmetry 
is explicitly broken if 6*2 7^ 0. The temporal component of the gauge field is diagonal in the flavor 
space, because the color and the flavor spaces are completely separated out in the present case. In 



the Polyakov gauge, L can be written in a diagonal form in the color space [|29|| 



where 0^ = ^3 + (ps/V^, (pb = -h + (ps/V^ and (pc = -{(pa + h) = -2(ps/V3. The Polyakov 
loop ^ is an exact order parameter of the spontaneous Z3 symmetry breaking in the pure gauge 
theory. Although the Z3 symmetry is not an exact one in the system with dynamical quarks, it 
may be a good indicator of the deconfinement phase transition. Therefore, we use <P to define the 
deconfinement phase transition. For simplicity, we assume below that mu = rrid = mo- 
We transform the quark field q to the new one q' with 

Qr = e'^QR, Ql = e-'^q'L, (6) 

in order to remove the 9 dependence of the determinant interaction. Under this U(1)a transforma- 
tion, the quark-antiquark condensates are also transformed as 

a = qq = cos {^^' + sin (0^'' 
V = ?«75g = ~ sm I - 1 a + cos I - Ir] , 
ai = qnq = cos I - I + sm I - I vr^, 

TTi = qira^q = - sin a- + cos 7i[, (7) 
where a' is defined by the same form as a but q is replaced by q'; this is the case also for other 
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condensates i]', a[ and vr-. The Lagrangian density is then rewritten with q' as 



C = q'it^,D''-mo+-imo--f,)q'~Ui<P[A],<P[A]\T) 

3 

+ G,J2[{q'T^qr + {q'^l,Taqr] 



a=0 



+ 8G2 [det (q'WL) + det [q'M] 



(8) 



where mo+ = mo cos (|) and mo- = mo sin (|). Making the mean field approximation and 
performing the path integral over the quark field, one can obtain the thermodynamic potential i? 
(per volume) for finite T and yu: 



(27r)3 



3{E+(p) + E_(p)} + - In [1 + 3^e-('^+ + 3^*6'^^^+ + e-^^^~^ 



+ 1 In [1 + 3<Pe-^^- + 3^*e-2^^- + e'^^^- ] + ^ In [1 + 3^*e-^''t + S^e-^/^^i + e-^'^^t] 



+ ^ In [1 + 3<P*e-^^~ + 3^e-2/3£;l ^ ^-3/3^1] + u + U, 
where i?^ = -E'+(p) ± /i and i?^ = (p) ± /i with 



(9) 
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p2 + C ± 2/D, 
+ A^^ + + 
A^M^ + P^N^ + 2APMN cos cp + A^P^ sin^ ip 
(MA + A^P)2 + (A X P)2 > 
mo+ - 2G+a = mo+ - 2{Gi + G'2)o-', 
mo- - 2G-r]' = mo- - 2{Gi - G2W, 
{-2G-a[, -2^-4, -2^-4), 
(-2G+7ri, -2G+7r2, -2G'+7r3), 
VA- A, P = VP P, A • P = cosy), 



/2 , /2n 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



In the right-hand side of dH), only the first term diverges. The term is then regularized by the 



1511 . we introduce c as Gi 



three-dimensional momentum cutoff A [|29, 3^. Following Ref. [|14l 
(1 — c)G and G2 = cG, where < c < 0.5 and G > 0. Hence, the NIL sector has four parameters 
of mo. A, G and c. We put mo = 5.5MeV. The parameters A and G are so chosen as to reproduce 
the pion decay constant f^^ = 93MeV and the pion mass m^ = 139MeV at vacuum. The remaining 



parameter c is a free parameter. Although the exact value of c is unknown, it is known from the 
analysis of the r]-r]' splitting in the three-flavor model that c ~ 0.2 is favorable pOll . The value 
c = 0.2 has been also used in Refs. [14, isj]. Therefore, we take c = 0.2 also here. 

The PNJL model is reduced to the NJL model in the limit of T = 0, since in Q the Polyakov- 
loop potential W tends to zero and the logarithmic terms approach NcO{ii — E^), where 0{x) = 1 
for X > and 0{x) = for x < 0, in the limit. Since P breaking at T = and 9 = n has already 



been studied in detail by the NJL model 113 
at finite T and 9 = n. 



15h . we here concentrate ourselves on P restoration 



The Polyakov potential U of Ref. 113 ih is fitted to LQCD data in the pure gauge theory at finite 



Ttl 



U = T^ 
a{T) 



a{T) 



+ h{T) ln(l - Q<P<P* + 4(^^ + <P*^) - 3{^<P 



2 

ao 



(19) 
(20) 



The parameters included in U are summarized in Table I. The Polyakov potential yields a first- 
order deconfinement phase transition at T = Tq in the pure gauge theory. The original value of Tq 
is 270 MeV evaluated by the pure gauge lattice QCD calculation. However, the PNJL model with 
this value of Tq yields a larger value of the transition temperature at zero chemical potential than 
the fuU LQCD simulation |53-3] 

predicts. Therefore, we re scale Tq to 212 MeV Il36ll so that the 
PNJL model can reproduce the critical temperature 173 MeV of the full LQCD simulation. 



ao 


ai 


02 


b3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential sector used in Ref. Bill . All parameters 
are dimensionless. 



The variables X = <P, <!?*, a, rci, r] and aj satisfy the stationary conditions, 

df2/dX = 0. (21) 

The solutions of the stationary conditions do not give the global minimum f] necessarily. There is 
a possibility that they yield a local minimum or even a maximum. We then have checked that the 
solutions yield the global minimum when the solutions X{T,9, jj) are inserted into (|9]). For 6* = 
and TT, the thermodynamic potential Q is invariant under P transformation, 

Tj -> -r/, Tla -TTa- (22) 
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Since the four-quark coupling constant G contains effects of gluons, G may depend on $ 



In fact, recent calculations [157 



5911 of the exact renormalization group equation (ERGE) [|6Qh 



suggest that the higher-order mixing interaction is induced by ERGE. It is highly expected that 
the functional form and the strength of the entanglement vertex G{<P) are determined in future by 
these theoretical approaches. In Ref. [42], we assumed the following form 



G{<P) = G[l - ai<P'P* - a2i<P'-^ + ^*^)] 



(23) 



by respecting the chiral symmetry, P symmetry, C symmetry [|l6l BSll and the extended Z3 sym- 
metry. This model is called the entanglement PNJL (EPNJL) model. The EPNJL model with the 
parameter set [42I, ai = ao = 0.2 and Tq = 190 MeV, can reproduce LQCD data at imagi- 



zero chemical potential [|53 



nary chemical potential |61-70] and real isospin chemical potenti al p71i1 as well as the results at 



55ll . Recently, it was shown in Ref. [|43ll that, also under the strong 



magnetic field, the EPNJL model yields results consistent with the LQCD data tlM- The EPNJL 
model with this parameter set is also applied to the present case with real /x and finite 9. 



III. NUMERICAL RESULTS 



The present PNJL model has eight condensates of quark-antiquark pair. However, a and vr 
vanish lll4l - ll6ll . since = and the isospin chemical potential is not considered here. We can 
then concentrate ourselves on a, rj and (p. We can also restrict 6* in a period < 6* < 27r without 
loss of generality. 

Figure [H shows T dependence of a and ^ at = and 6* = 0; note that 77 = in this case. 
In both the PNJL and EPNJL models, a rapidly but continuously decreases (increases) as T 
increases. Therefore, the chiral restoration and the deconfinement transition are crossover. The 
crossover transitions occur more rapidly in the EPNJL model than in the PNJL model. 

Figure |2] presents T dependence of a, rj and <P at /i = and 9 = n; note that this figure plots 
condensates a and i] before the transformation (|6]). For both the PNJL and EPNJL models, the eta 
condensation {i] ^ 0) occurs at low T and hence P symmetry is spontaneously broken there. At 
high T, in contrast, i] = and hence P symmetry is restored. The critical temperature Tp of P 
restoration is 202 (170) MeV in the PNJL (EPNJL) model. 



The order of P restoration was reported to be of second order in the NIL model [|l4l 



Jj but 



of first order in the linear sigma model [|73il . In the EPNJL (PNJL) model, r] is discontinuous 
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Fig. 1 : T dependence of the chiral condensate a (solid line) and the Polyakov loop <P (dashed line) at ^ = 
and 9 = 0. Panel (a) represents results with PNJL model, while panel (b) does results of the EPNJL model. 
The chiral condensate is normalized by the value uq at T = 9 = 0. 

in its zeroth (first) order. Thus, the PNJL model supports the second-order transition, but the 
EPNJL model does the first-order. The zeroth-order discontinuity (gap) of 7] in the EPNJL model 
is propagated to other quantities a and ^ as the zeroth-order discontinuity (gap), according to the 



discontinuity theorem by Barducci, Casalbuoni, Pettini and Gatto [17411 . This really takes place 
in Fig. |2l although the discontinuity (gap) is appreciable for ^ but very tiny for a. The first- 
order discontinuity (cusp) of r] in the PNJL model is also propagated to a and <1> as the first-order 
discontinuity (cusp) [37]. As mentioned above, the EPNJL model is more consistent with the 
LQCD data than the PNJL model. This means that the order of P restoration may be weak first 
order. 

T dependence of cr, r] and <P is shown also for /i = SOOMeV and = tt in Fig. [3l In the case of 
large fi, P restoration is of first order in both the PNJL and EPNJL models. As jj, increase with 
fixed at TT, thus, P restoration changes from the second order to the first order in the PNJL model, 
while it is always of first-order in the EPNJL model. This means that there is a tricritical point 
(TCP) in the PNJL model; note that the TCP is a point where the first- and second-order transitions 
meet each other. 

Figure |4] shows the phase diagram of P restoration in the ji-T plane for the case of = tt. In the 
PNJL model of panel (a), point A at (/x, T) = (209 [MeV], 165 [MeV]) is a TCP of P restoration, 
while there is no TCP in the EPNJL model of panel (b). 

Figure |5] represents the phase diagram of the chiral transition in the n-9-T space. These dia- 
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Fig. 2: T dependence of the chiral condensate a (solid line), the eta condensate i] (dotted line) and the 
Polyakov loop (p (dashed line) at jj. = and 6 = n. Panel(a) stands for results of the PNJL model, while 
panel (b) does results of the EPNJL model. 




Fig. 3: T dependence of the chiral condensate a (solid line), the eta condensate r] (dotted line) and the 
Polyakov loop (p (dashed line) at fi = SOOMeV and 9 = tt. Panel(a) stands for results of the PNJL model, 
while panel (b) does results of the EPNJL model. 



grams are mirror symmetric with respect to the fi-T plane at 6' = tt. At finite 9, a/ct^ + 77^ is the 



order parameter of the chiral transition rather than a itself [11611 . In the PNJL model of panel (a), 
point A in the fi-T plane at 6* = tt is a TCP of P restoration and a critical endpoint (CEP) of chiral 
restoration at which the first-order (solid) line is connected to the crossover (dotted) line. Point 



75|1 . The second- order 



C in the fi-T plane at 6* = is another CEP of the chiral transition [|46 . 
(dashed) line from C to A is a trajectory of CEP with respect to increasing 9 from to tt. Thus, the 
CEP (point C) at 6* = is a remnant of the TCP (point A) of P restoration at 6* = tt. In the EPNJL 
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Fig. 4: Phase diagram of P restoration at = vr. Panel(a) stands for results of the PNJL model, while panel 
(b) does results of the EPNJL model. The solid and dashed lines represent the P restoration of first-order 
and second-order, respectively. Point A in panel (a) is a TCP 

model of panel (b), no TCP and then no CEP appears in the yU-T plane at9 = n. The second-order 
(dashed) line starting from point C never reaches the /i-T plane at 6^ = tt. For both the PNJL and 
EPNJL models, the location of CEP in the /i-T plane moves to higher T and lower ^as 9 increases 
from to TT. Particularly in the EPNJL model, the chiral transition is always first order in the yU-T 
plane at = tt. 




Fig. 5: Phase diagram of the chiral restoration in the — 9 — T space. The solid, dashed and dotted lines 
represent the chiral transitions of first-order, second-order and crossover, respectively. Point A is a CEP of 
chiral restoration and a TCP of P restoration, while point C is a CEP of chiral restoration. Panel(a) stands 
for results of the PNJL model, while panel (b) does results of the EPNJL model. 

Figure |6] shows the projection of the second-order chiral-transition line in the ix-O-T space on 
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the fi-O plane. The solid (dashed) line stands for the projected line in the EPNJL (PNJL) model. 
The first-order transition region exists on the right-hand side of the line, while the left-hand side 
corresponds to the chiral crossover region. The first-order transition region is much wider in the 
EPNJL model than in the PNJL model. In the EPNJL model, eventually, the chiral transition 
becomes first order even at = when is large. 
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Fig. 6: Projection of the second-order chiral-transition line in the fi-9-T space on the plane. The solid 
(dashed) line corresponds to the line in the EPNJL (PNJL) model. 

QCD has the anomalous Ward identities among the chiral and eta condensates, a and r], the 
gluon condensate {-^;^^°'^'^'^ F^yF^p) and the topological susceptibility Xt MM'- 



'de 



'647r 



> av a pi 



AT, 



= -Xt = -Trt^^oo- + 0{ml 



(24) 
(25) 



The identities are useful for checking the self-consistency of the proposed model [|l2n . The gluon 
condensate does not appear explicitly in the PNJL model, but the model satisfies 



(26) 



The PNJL model also satisfies the second Ward identity (l25l) as shown below. In the PNJL model, 
the left-hand side of (l25l) is rewritten into 



with the inverse curvature matrix 



fixed 



C5' 



E 



d'^Q ( d^Q d^Q 



ded(pi \d(f) 



(27) 



d(f)id(f)j 
11 



(28) 



for the parameters cj) = (a', i]', <P, <P*). Note that C~j^ is the susceptibility Xij of order parameters 
(pi and (pj and the term including Xij is of order 0{ml). Equation (ITTI) tums out to be the second 
Ward identity (|25l) after simple algebraic calculations. For 9 = 0, Eq. (|27l) is further rewritten into 



06^ 4 



-nioa + - 



12° 



4(2G'_f2^^ + l) 

_ 92^20 



0; 



(29) 



in the virtue of the 6'-reflection symmetry, where f^j^j^ = for the thermodynamic potential 
12° with no the meson potential U and no the Polyakov potential U. 
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Fig. 7: T dependence of the topological susceptibility xt (solid line) and the chiral condensate o (dashed 
line) at (a) = 0, (b) 7r/2, and (c) tt. These are results of the PNJL model at = 0. 

Figure |7] shows the topological susceptibility xt and the quantity mocr/4 as a function of T at 
6* = 0, 7r/2 and tt, where the case of = is considered. As shown in panels (a) and (b), the 
topological susceptibility almost agrees with mocr/4 for < 7r/2. The small deviation between the 
two quantities shows that the corrections of order 0{vn?^ to the approximate identity Xt = ]^"^oC" 
are small. In panel (c) for = 71, the both quantities are still close to each other except for the 
critical temperature T ^ 200 MeV of the second-order P transition. Near the critical temperature, 
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the susceptibility Xij of the order parameters becomes large and hence the corrections of order 
0{ml) coming from the second term of the right-hand side of (|T71 become significant. 

IV. SUMMARY 

In summary, we have investigated theta- vacuum effects on QCD phase diagram, using the PNJL 
and EPNJL models. For the both models, the chiral transition becomes strong as 9 increases. 
Particularly in the EPNJL model that is more reliable than the PNJL model, it becomes first order 
even at /i = when 9 is large. This is an important result. If the chiral transition becomes first 
order at zero yU, it will change the scenario of cosmological evolution. For example, the first- 
order transition allows us to think the inhomogeneous Big-Bang nucleosynthesis model or a new 
scenario of baryogenesis. Our analyses are based on the two-flavor PNJL (EPNJL) model and 
effects of the strange quark are then neglected. It is very interesting to study the effects by using 
the three-flavor PNJL model [|40|]. 
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